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Abstract 

We found the eigenvalues of the transfer matrices for the 1-D Hubbard model 
and for the coupled XY model with twisted boundary condition by using the 
analytic Bethe Ansatz method. Under a particular condition the two models 
have the same Bethe Ansatz equations. We have also proved that the periodic 
1-D Hubbard model is exactly equal to the coupled XY model with nontrivial 
twisted boundary condition at the level of hamiltonians and transfer matrices. 
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1 Introduction 



The 1-D Hubbard model (the Hubbard model) is one of the significant exactly solvable 
models in condensed matter physics. Lieb and Wu [jT| succeeded in diagonalizing the 
hamiltonian in the frame of the coordinate Bethe Ansatz. However, the integrability 
of the Hubbard model was recently set up by Shastry []2], |J, Olmedilla, Wadati and 
Akutsu 0, |5|, |6| from the viewpoint of the Quantum Inverse Scattering Method (QISM). 
The Yang-Baxter equation for the related R matrix was proved in 0. In [§, H|, Shastry 
found the Yang-Baxter relation and constructed the transfer matrix of the coupled XY 
model which is equal to the Hubbard model with the help of the Jordan- Wigner trans- 
formation. The commutative family with one free parameter ensures the integrability 
of the system. In a different approach, Olmedilla, Wadati and Akutsu 0, ||, starting 
from the super L-operator of the Hubbard model, solved the super- Yang-Baxter (SYB) 
relation and found the invertible R matrix which is the same as that given by Shastry 
P| up to a scalar function. In both cases || U, the hamiltonian can be derived from 
the transfer matrix under the periodic boundary condition. 

The base of the Quantum Inverse Scattering Method is the Yang-Baxter relation. 
The later has closely related to the Quantum group and the Yangian. This technique 
provides a systematic method to deal with the integrable 1-D quantum systems and 
2-D solvable statistical mechanical models. It is known that most of integrable systems 
can be handled in the frame of the Algebraic Bethe Ansatz ( or Analytic Bethe Ansatz 
) method. But, up to now, there is no report about the diagonalization of the transfer 
matrix for the Hubbard model in the QISM approach. So, it is important to find 
the solution of the Hubbard model by using QISM. In ||, Shastry conjectured the 
eigenvalue of the transfer matrix of the coupled XY model based upon the coordinate 
Bethe Ansatz method. Comparing the Bethe Ansatz equations given by Lieb and Wu 

and by Shastry H, an extra factor appears and it shows some difference between 
the two models. 

The motivation of this paper is to find the eigenvalues of the transfer matrices 
related to the Hubbard model and the coupled XY model with twisted boundary con- 
dition by a version of analytic Bethe Ansatz method (ABA). We want to give some 
explanation about the difference between the coupled XY and the Hubbard models. 
The extra factor originates from the boundary condition. In fact the Jordan- Wigner 
transformation does not keep the invariance of the boundary condition. We also show 
that the Hubbard model with periodic boundary condition is exactly equal to the 
coupled XY model with a special boundary condition. 

The organization of this paper is to recall the Yang-Baxter relation for the Hubbard 
model and the parametrization of the R matrix in section 2. In section 3, we investigate 
the model at special cases by using the algebraic Bethe Ansatz method. Because there 
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are two kinds of creation operators with spin-up ( or spin-down ) , such as T21 and T43 
(T 2 i and T 43 ), the general multi-particle states become very complicated. However, the 
special solution gives an insight about the general structure of the eigenvalue. Section 
4 will attributes to an analytic Bethe Ansatz method to the Hubbard model. In our 
approach, the analytic property of the eigenvalue of the transfer matrix, together with 
the asymptotic behaviour, determines almost all unknown functions. The discussion 
is different from the standard analytic Bethe Ansatz method since there is not the 



crossing symmetry, which has played an important role || [10], [11], |l2| . In section 5, we 
will apply the analytic Bethe Ansatz method developed in section 4 to the coupled XY 
model with twisted boundary condition. Under a special choice, it recovers the results 
of the Hubbard model in the level of Bethe Ansatz equations and the eigenvalue. In 
section 6 we show that the Hubbard model with periodic boundary condition is equal 
to the coupled XY model with special boundary condition by using the Jordan- Wigner 
transformation and some discussions are given in section 7. 



2 The Hubbard model and the Super- Yang-Baxter 
relation 

In this section we recall the definition of the SYB relation for the Hubbard model and 
some useful functional relations. We follow the notation in ||. The hamiltonian of the 
Hubbard model is 

L L 

H Hu = - Yl ( a m+i,s a m,s + a^ s a m+1 ^ +Uj2{n m ^ -l/2)(n ml - 1/2) (1) 

j=l,s=1,l m=l 

where a+ s (a mjS ) stands for the m-th site electron creation ( annihilation ) operator 
with spin s. The Super- Yang-Baxter relation is || 

K(n, v) [£ m (/i) ® s C m {v)\ = [C{y) m ® s C m (n)} K(p, v) (2) 

The super tensor product is defined by 

[A ® s B]f b = A}B^(-1)^ +Pij ^ a \ (3) 

where p(l) = p(4) = —p{2) = — p(3) = 1 are the parities. For the Hubbard model the 
L-operator takes the form 

£(Ai) = 
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ml 



-h(n) 



ia m1 a ml e hM \ 



^/m| a m | fm\9ml& ^ 

a ml fml 

ia m\ a ml^ Kil) a m\9ml Wm^l ~9m\9mlZ 

where \i is the spectrum parameter, and the functions f ms and g ms are 
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9ms 



= Wi(n) - w 3 (fi) - {w 4 (ii) - w 3 (fi) - i[wi(n)) + w 3 (/j,)]}n ri 



w 4 (n) +w 3 (n) - {w 4 (n) 
w A (fx)+w 3 (fi) ■ 

sinh(2/i(/i)) : 



w 3 (fi) - i[w 4 (ii) - w 3 (n)\}n ri 
sin(yU + 7r/4) 

sin(yU — 7r/4) 
--cos(2/x) 



The graded 7?. matrix reads 
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where all functions are defined by 

p 4 (/i,z/ 
p 9 (p,u 

Pe(p,u 
p 7 (p,u 
Pz{p,v 



e l a(p)a(u) + e~ l -f(p)-f(u), 

e l l{p)l{v) + e~ l a(p)a(u), 

-e l a(p)-f(v) + e~ l -f(p)a(u), 

e l l{p)a{v) - e~ l a(p)-f(u), 

e l a{p)a{v) — e~ l ~f{p)~f{v) 
a 2 {p) - 7 2 (z/) 

-e lr )(p)^(v) + e~ l a(p)a(v) 
a 2 {p) - 7 2 (z/) 

e~ h [e l a(p)^(p) - e~ l a(u)^(u)] 

a 2 {p) -7 2 (^) 
Pi(p, v) - p 3 (p, v), 



(7) 



= pi(n,v) -p 5 (p,u), 

I = h(p) — h(u), 
h = h(p) + h(u). 

Due to the Super- Yang-Baxter relation, one can define the monodromy matrix 

T(p) = C L (p)---C 2 (p)C 1 (p) (8) 

which still satisfies the Super- Yang-Baxter relation 

TZ(p, v) [T(p) ® s T{y)\ = [7» ® s T{p)\ TZ(p, u) (9) 

The Super- Yang-Baxter relation leads to the existence of the commutative family of 
the transfer matrices tn(p) = strT(p) with infinitely many different value of p. So, 
the Hubbard model is integrable. The infinite number of conserved quantities can be 
derived from the transfer matrix The derivation of log[t#(/x)] at p = tt/4 gives 

the hamiltonian of the Hubbard model with the periodic boundary condition. Before 
ending this section, we list some useful functional relations which will be used in the 
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following sections, 



1 


= PlP4+P9PlO, 




2 


= PlP5 + P3P4, 
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= PzPh - (Pe) 2 , 




PlO 


= pQ(e h a(n)a(p) + e 


- h i(ph(")) 


P9 


= P6(e~ h a([i)a(v) + 


e h l(p)l(v)) 



(10) 



3 An algebraic analysis on the eigenvalue 

In this section, we will discuss the solution of the Hubbard model in some special cases 
by using the algebraic Bethe Ansatz method. For the general multiparticle states, it 
becomes very difficult even if there are two electrons with opposite spins. 

Taking the special elements of the Super- Yang-Baxter relation, one can obtain the 
following useful relations 
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(13) 



T 11 ( /U )T 43 (z/) = 



nfl(fi,p) 

KfM v ) 

n%(^p 



T 21 (n)T M (p) 
T A2 (fi)T 23 (p) 



TZfl(ii,iy 
TZfKn, p 

nfM v 



T 43 (p)T 11 ( f i) + 



n?M») 



V 41 (n p) 



(14) 



KfM») 

For the Hubbard model the Hilbert space consists of four states: double occupied 
state | ][), spin-up state | ]), spin-down state | [) and unoccupied state |0). We denote 
them by |1), |2), |3) and |4) respectively. It is convenient to introduce the reference 
state 

\vac) = |4)i ®i • • • <g) s |4) L (15) 

Using the explicit expression of the L-operator, one can find that the monodromy 
matrix acting on the reference state takes the form 

/ Ai(/i) \ 

T 21 (/i) AM 
T 31 (/i) AM 



T(/i)\vac) 



\vac) 



(16) 



where 



\ %M T A2 {y) %M M{n) ) 
AM) = i-aMe hM ] L 
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AM = AM) = [a(/i) 7 (^)e-^] L 

AM = [- 7 2 (/i)e ft ^] L 

Unlike the case where the nested Bethe Ansatz method is applicable, the operator 
is related to operator 743, creating an electron with spin-down. Generally, the 
relation is very complicated. When they act on the reference state, however, the 
situation becomes simple. An algebraic calculation shows 

where the quantity with tilde is defined by 

e -2fc(A)^) =e 2/M^) (lg) 

7(A) i{p) 

Further, one can show 

72i(a0 743(^1) ' • -TiziUn^vac) oc T^/^T^/ii) • • ■T 43 (fj ln )\vac). 

It is worthy to notice that this relation is valid only for the reference state; when acting 
on the other states, it will be invalid. There is a similar relation between %i and T 42 . 
Due to this relation, we can construct the special states with all spin-up (spin-down) 
with 71 3 (7^2) 

Using the commutative relations, we can find 



where 



tH{n)Y&N) = A(/i) I ^> N ) + unwanted terms. 



m = am n fl r 2 A - m*) n . ^ 

j=1 ipM'H j=i wMm 

-am n ^10(^^0 

uwn -7^7) (Mi 

jJiPM^j) - pMPj) 



The vanishing unwanted terms gives the Bethe Ansatz equation 



[-e 2 w^p(] L = 1 (20) 
l\Pj) 
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Thus, the states \^n) are really the eigenstates of the transfer matrix t(p) if the 
spectrum parameters are appropriately chosen to satisfy the Bethe Ansatz equation 
(20). 

The general states with iV spin-down and N — M spin-up are given by sums of 
products of the combination of T 2 i, %i, and T^,j = 1,2,3 acting on the reference 
state. Let us consider a special case N — M — M — 1. The general form is 

= {/i (^1,^2)^42(^1)^43(^2) + Mvi, fa)%aM%2M 

+ /3(A*1, ^2)^21(^1)^31(^2) + /4^1, ^2)^3l(/^l)^2l(/^2) 

+ h{pu P2)%i{pi) + /e(A*i, A*2)^4i(a*2)} \vac) (21) 

We can determine the all coefficients fj by requiring the r.h.s. of equation (21) to 
be the eigenvector of the transfer matrix. The general solution is very complicated. 
Fortunately, we can show by an explicit calculation that is the eigenstate if 

/3 — fi — with appropriate / 5 and f$. Using the Super- Yang-Baxter relation, we 
find the eigenvalue for the case when f\ = —f 2 
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-ipio(p, Hj) 



-Mp) 



,=i itoiHj, p) jA Pi (A*, Pj) ~ Pj) 

p 4 (/i,/il)pi (/i,/i 2 ) 



\p 9 (p,/il)[pi(/i,/i 2 ) - p 3 (p,ll2)} 

Pi{PiP2)Pw{PiPi) 



Pq(a*> A*2)|>i(a*> A*i) -Ps(^,/^i)] 



Pw{p,Pi 
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Pl(p,pl) - p 3 (p,pl) 

Plo(/i,/i2) 
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P${P,Pl) 

Pa{p,P2, 



_Pi(p, P2) - Ps(p, P2) pg(p, P2)_ 
The vanishing unwanted terms gives the Bethe Ansatz equation 

~ —a(p^ 



2/i(/i, 



)] L = 1 



(22) 



(23) 



For the case /1 = / 2 the eigenvalue and the Bethe Ansatz equations are given by 
the equations (19) and (20), respectively. It is worthy to point out that these states 
considered here are not complete. For example, consider the case f\ — f2 — 0, one can 
get similar results. 
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4 Analytic Bethe Ansatz for the Hubbard model 



In the last section, we applied the algebraic Bethe Ansatz method to some eigenstates 
of the Hubbard model. For general states, the straightforward calculation becomes 
very complicated. In this section, however, we want to discuss an analytic Bethe 
Ansatz method to the same problem based on the hints given by the above results. 



We should generalize the standard ABA [BL [TH, O, [L2| in which the crossing symmetry 



and asymptotic behaviour play a key role. We can not apply the same argument to 
the eigenvalue of the Hubbard model in which there is no such crossing symmetry for 
the R matrix. 

Let us first investigate the special solutions and show how to generalize it from 
the viewpoint of ABA. One may understand that the analytic property of function 
A(/i) leads to the Bethe Ansatz equation (20). It is clear that Pq(p, f-ij) = is the 
simple pole of A(^) (equation (19)). In order to keep the analytic property of the 
eigenvalue, the residue at such pole must be zero. The Bethe Ansatz equation is 
nothing but the condition of vanishing residue. Similarly the vanishing residue at the 
pole pi(/i, /ij) = ps(p, Pj) gives the same Bethe Ansatz equation. This property can 
be generalized to all kinds of states with different particles. The eigenvalue function 
should be analytic and has only superficial simple poles. The vanishing residues at 
such poles will give the Bethe Ansatz equations. 

Let us discuss the general eigenvalue of the Hubbard model. The special solutions 
(19) and (22) together with some standard knowledge of algebraic Bethe Ansatz contain 
enough message about the general one. It consists of four terms which are proportional 
to Aj(fi),j = 1, 2, 3, 4 respectively. The terms involving A\ and A 4 are dependent only 
on the total number of electrons. The other terms depend on both the total number 
N of electrons and the number M of spin-up electrons as shown in the last section. 
Thus, the general eigenvalue should be 

m = ma n - am ft ~ p f^ ^ n ^ \ m ) 

N _■ ( \ M 

-am n i \ n nit*, 

+ A M fl -'»^>f (24) 

j=l PlKPiPj) ~ P3[P,Pj) 

where #2 arid g% are undetermined functions. The pj and A m are free parameters. N 
is the total number of electrons, M the number of spin-up electrons. 

Now, we show how the analytic property of the eigenvalue restricts the unknown 
functions. First, A(/i) has two sets of poles related to parameters fij. One (case A) is 
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controlled by the null denominator of the first two terms of equation (24). The another 
(case B) is from the last two terms. For case A the position of poles is determined by 



a(p) 
70) 



(25) 



Due to the 27r-period of h(p), we can get /x = fij in the region < /i j < ir. At these 
poles, the functions pi(/i, Pj) — ps(p, fij) and pw{p, Pj) also vanish, But, the ratio is 
finite. So, the singularity at these poles is dominated by the terms in the first line. 
The null residue requires 



3 2/t( W ) 



a(nj 



L M 



(26) 



m=l 



For case B, the position of poles satisfies 
= 



(27) 



The first part is equal to pg(p,pj) = 0. It is not a pole due to pw(p, Pj) = at this 
point. The real pole locates in jlj 



a(flj 



a(pj 



(28) 



The vanishing residue at /i = pj gives 



_ e 2h{fij) 



L N 

= rifi , 2(/ij,A m ) 

3=1 



(29) 



In order to keep the analytic property of the eigenvalue, Equations ( p6|) and (|28|) must 
be satisfied simultaneously, which leads to the following functional relation 



Ai 



M 



m=l m=l 

We find that it is convenient to introduce the new variables k and kj 

Jk _ 2fe(^) Ct (/ i ) 



(30) 



— e 



7(/^) 



(31) 
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In terms of the new variables the relation between the quantities with and without 
tilde is very simple 

sm(kj) = sin(kj) + i— (32) 

At this stage, we need to know some properties of the undetermined functions 
which hides already in the special solution with N = 2, M = 1. In terms of kj, equation 
(f22|) changes into 

A(ji) = {-a 2 (/i)e h ^) L A(fc) 

A = ^ 2icos((fc + fc J -)/2)co S ((fc + fc J -)/2) 

i sin (A;) — i sin(kj) 



N 

n 



lkL SL 2i cos((k + kj)/2) cos((k + kj)/2) 



. , % sin (A;) — z sin(fcj) 

i2 sin (Ac) — i sin(iTi) — i sin(&2) — U/2 
i2 sin (A;) — z sin(/ci) — i sin(/c 2 ) + U/2 



2 



2% cos{{k + kj)/2) cos((A: + kj)/2) 



- e~ ikL TT 

f-Jl i Bjn(Jfe) - i sin(fcj) + [7/2 

z'2 sin(A;) - % sin^i) - z sin(fc 2 ) + 3E7/2 
z'2 sin(A;) — z sin(A; 1 ) — z sin^) + U/2 

_ p -i(k+k)L A 2icos((fc + fc/)/2)cos((fc + fc/)/2) 

^ i sin(ifc) - i sin(%) + C//2 1 J 

It is clear that the eigenvalue has another simple pole 2 sin (A;) = sin(fci) + sin(A;2) — U /2. 
The Bethe Ansatz equation ensures the analytic property of the eigenvalue. This 
strongly suggests that the undetermined functions have the simple poles of the form 
sin(fc) = const. On the other hand, the analytic property of the eigenvalue requires 
that g 2 and g% must have same poles. Therefore, the general form of ^(a*, A) and 
9a{p, A) is 

P 2 (M) ( ,s P 3 (M) 

z sin(fc) — A + (7/4 z sin(fc) — A + (7/4 

where the function P 2 (/c, A) and Ps(k,X) are integral functions. The general form is 
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P r (k,X) = J2n=o a n(k)W n - Substituting it into equation (30), we have found 

a 2 (k) = [ism(k)-U/A]al(k) 

a 3 (k) = [ism(k) + 3U/4](al(k) 

(35) 

a 2 n (k) = al(k) = 0,n>2 

af(k) =[aj(k)}^ 

where k is defined by sin(A;) = sin(A;) — iU/2. Moreover, the function a\(k) is analytic 
and has no zero in the complex plain. It will be fixed by the asymptotic behaviour 
of the transfer matrix. Let us assume U < and fi — > —ioo, then the eigenvalue 
approaches to 

A(/i) — e^{[aj(oo)] M + [a?(oo)]- M }. (36) 

Comparing this with the asymptotic behaviour of t(fi) — > e 3Lo °, we get af(oo) is a 
no- zero constant. Based upon the knowledge of analysis such as Liouville's theorem on 
integral functions, we arrive at a\{k) being a no-zero constant. A special form of A(/x) 
under N = M = 1 fixes this constant to be unit. Finally, we arrive at the final results 

m = (-« 2 (/i)e^) L A(k) 

= ^ 2tcos((fc + fc 3 -)/2)cos((fc + fc J -)/2) 
jj[ isin(k) — isin(kj) 

_ e - lkL fr 2tcos((fc + fc J -)/2)cos((fc + fc J -)/2) 
j-Ji i sin (A;) — % sin(A; J ) 

^ ism(k)-\ m -U/A 
X ii 1 isin(A;)-A m + C//4 



N 



-ikL 



2% cos((A; + kj)/2) cos((A; + k j )/2) 



-e """TT 

jJi i s in(A;) - i sin(^) + U/2 

™ ism(k)-X m + 3U/A 

x n 



m- 



^ i sin (A;) — A m + C//4 



_ p -i( k+ ~k)L A 2^cos((A: + A: J )/2)cos((fc + A: J )/2) 

^ i sin (A;) - i sin(A;.,) +U/2 { ' 

The parameters satisfy the following Bethe Ansatz equations 



e 



ikiL 
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n 



m- 



i sm.{kj) — \ m — U/A 
\ i sin(/c J ) — \ m + U/A 
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__ -ft Xr-X m - U/2 " zsin(fc J )-A r + ^/4 

ii A r - A m + U/2 " Jti i sinCAy-) - ^ - Z7/4 1 j 

Differentiating log(A(/x)) at // = 7r/4 will give the energy of the Hubbard model, which 
coincides with the one given in |J 

£=- ^-£ cos (%) (39) 

* z i=i 

We have checked that equations ([37]) and (|38| ) with AT = 2, M = 1 coincide with the 
result obtained in section 3. 



5 ABA for the coupled XY model 



In this section we investigate the eigenvalue of the transfer matrix of the coupled XY 
model with twisted boundary condition. 

The L-operator related to the coupled XY model is [|]. [| : 



(40) 



where cx^ and are two independent Pauli matrices located in m-th site. The operator 
•p^ and q ± read 

(41) 
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In references @, , it was shown that this L-operator satisfies the Yang-Baxter equation 
R(jjl, v)L m (n) <g> L m {y) = L m (u) <g> L m (ji)R(ji, u) (42) 
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The R matrix is 
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(43) 



This Yang-Baxter relation ensures the monodromy matrix T(/x) = Li(fi) ® • • • <g> Za(p) 
satisfying the Yang-Baxter relation. In order to simplify our calculation, we choose the 
ferromagnetic state ( all spin-down states) as the reference state. From the explicit 
expression of the L-operator, we have 



T(n)\vac) 



( ^i(p) \ 

T 2 i(p) A 2 (n) 

T 31 (/i) AM 

^ T 4 i(/i) T 42 (/i) T 43 (/i) Ai(/i) J 



\vac) 



(44) 



where 



AM = AM = K/i) 7 (p)e-^] L 

AM = [7 2 (p)e^ } ] L (45) 

Similarly one can define the transfer matrix t(/i) = stT(fi) and find the eigenvalue 
of it, which will related to the periodic boundary condition. The eigenvalue of the 
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diagonal-of-diagonal transfer matrix of this model with periodic condition was found 
by Bariev [|13[ in terms of the coordinate Bethe Ansatz method. 

In order to consider twisted boundary condition, we introduce the generalized trans- 
fer matrix 

= T u ( M )a^ + T 22 ( M )a^ + TMa^ + T u (fj.)a ifk (46) 

where 



Jf3 2 



i(3 4 



o a a N a +a r N T +a 



^CcrNa+CrNr+CO 



(47) 



where a s and c s are free parameters, N a (N T ) the total number of cr-spin (r-spin). Now, 
we want to find the eigenvalue of t 9 (n) by means of the analytic Bethe Ansatz method. 
First, by using the algebraic Bethe Ansatz method, we find the eigenvalue of the states 
with N r-spin (or a-spin) flipping from the reference state. After a long but direct 
calculation, we arrive at 



Ajv(aO 



ee asN+ao {e hM a 2 (n)] L ll 



N 



Pl{Pj,P) 



r- 



IP^PjiP) 



+e' e c ^ +C0 [e- h ^a(/i)7(/i)] L fl ^ H) 



=1 Pv{P,Pj) 



N 

f e-^ N - C0 \e- h ^a(^^)] L ll 



Pw(p,Pj) 



(4£ 



+e e 



-a 3 N~a \ e h(n) 2 



N 



[e h{ »h\p)] L n 



iPi(Pj,p) - Ps(P,Pj) 
Pi(p,Pj) 



3=1 PsiPiPj) ~Pl(P,Pj) 

where s = a (s = r) for all a (r) spin-up states, the parameters /ij are determined by 

^2h{m) a \P'h 



' _]\N+1 e (c s -a s )N+co-ao 



(49) 



i(Pj)a 

^From this expression and the similar arguments in the last section, we can write the 
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general form of the eigenvalue 

A(/l) = e a +a„M+a T (N-M){ e h Ma 2^L -Q Pj^Wfj 

j=i ftlftiW 

+eC0+C(T M +CT (iV-M) [e -MM) a ( /i ) 7 ( /i )]i 



i=1 M/^jJ m=l 
+e -co-c CT M- Cr (iV-M) [e -MM) Q; ( // ) 7 ( // )]i 



(50) 



Second, consider the singularity of the A(/x) at the poles related to the parameters 
/ij. As done in the Hubbard model, the null residue condition requires the following 
relation 

M M 

n 9^Am) = n ft^oM™) (51) 

m=l m=l 

This equation is same as that in the Hubbard model (equation (34)). So, we can use 
the results of the Hubbard model. 



i sm(k) - A - U/A 
92[ ^ X) ~ C *sin(A;)-A + [//4 



A) = 



1 i sin(fc) - A - U/A 
cism(k) - A + C//4 



(52) 



here we have used the same definition of k as one in the Hubbard model. One should 
note that in this case, the constant c in the above equation not being 1 as in the 
Hubbard model. Taking N = M = 1 in equation(50) and comparing with equation 
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(48), one can get c = — 1. Thus, we obtain the final results 
A(jj) = [e h ^a 2 (fi)] L A(k) 

A r k \ = ( l\N c a +a a M+aAN-M) A 2 COS (( fc + k j)/ 2 ) COS ((k + kj)/2) 

j=i i sin(k) — i sm(kj) 

+eC0+CaM +c T (N-M) t~t 2 cos((fc + kj)/2) cos((fc + kj)/2) 
jJi i sin(k) — i sin(kj) 



i i \L+m —ikL fr ism(k) - X m -U/4 

' ^ 1 tsm(k)-X m + U/4 

A' 



+e _ C0 _ C(T M_ CT{A r_ M) ^ 2cQ S ((fc + fc j )/2) COs((fc + fcj)/2) 

^ i sin(fc) - z sin(%) + U/2 

y, I -i \L+M+N —ikL TT * sin(fc) - A m + 3^/4 

1 J ^ isin(k)-X m + U/4 

+(r a . aaM -a T {N-M) e -i(k + k)L -Q 2 COs((fc + kj)/2) CO S ((fc + fcj)/2) 



2 sin (A;) — i sin + U/2 



The Bethe Ansatz equations are 

^_^^A/+A r +l+-L e «fc J L _ e c -a +(c {T ~a a )M+(c T -a r )(N-M) 



x g isin^-A^-^ 



N 



ism(kj) - X r + U/4 



^isin^) - A r -[//4 ~ 1 



^ isin(kj) — X m + U/4 



N+l e 2(c +c a M+c T (N-M)) 



(53) 



(54) 



x n ^ r ~^ m ~^^ 

m=l ^ r ~ ^ m J rU/2 

These are the exact solution of the coupled XY model with twisted boundary condition. 
It is of interest that it recovers the results of the Hubbard model when a a = a T = 
c a = —c T = — in/ 2 and a = iLn,co = in. This gives a precise relation between the 
two models and make it clear why the extra factor appears among the Bethe Ansatz 
equations by Lieb and Wu |I| and by Shastry ||. When a s = c s = 0, s = a, r, 0, they 
reduced into the periodic case. The correspondence between our notation and that in 
H is 2ism(kj) = zj 1 — Zj. 
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6 Twisted boundary condition 

In this section, we will discuss the boundary condition related to our generalized trans- 
fer matrix (46) and prove the equality of the periodic Hubbard model to the twisted 
coupled XY model in terms of the transfer matrices and the hamiltonians. 

First, we derive the hamiltonian related to the t 9 (46) by using the standard method. 
After a straightforward calculation, we arrive at 

L-l rj L 

" = L^i \ a m+l a m + a m a m+\ T m+l T rn + T m T m+l) + ~T L^i <J rn T m 
m=l m=l 

+ exp{-e[a + c + (a a + c a )N a + (a r + c T )N T }}aj^a^ 

+ exp{e[a + c + (a CT + c a )N a + (a T + c T )N T ]}a N a+ ( 55 ) 

+ exp{-e'[a - c + (a a - c a )N a + (a T - c t )N t ]}t^t{ 

+ exp{e'[a - c + (a CT - c c )N a + (a T - c T )N T ]}r N r^ 

where 

{1 rspin up ,_J 1 aspm up , , 

— 1 rspin down ' | — 1 aspin down 

This means that the hamiltonian (55) gives the coupled XY model discussed in the 
last section with the twisted boundary 

<7± +1 = e" e(aCT+CCT) exp{±e[a + c + (a, + c a )A^ CT + (a T + c r )A^ T ]}(7 1 ± 

(57) 

= e- e ' {aT - CT) exp{±e[ao-Co + ( y a a -c a )N a + ( y a T -c T )N T ]}T^ 

Comparing the Bethe Ansatz equations (37), (38), (53) and (54), we find that they are 
same if the free parameters are fixed 

a a = a T = c a = — c T = — — , ao = iLn, cq = iir (58) 

Let us prove the connection (58) by using the Jordan- Wigner transformation on 
transfer matrix. The Jordan- Wigner transformation on operators is defined by 

' " + N = * ( t ) = ( 1" ) ( t ) (59) 



® m 



T m \ T/2 [ a ml 

t~, ~ ml \ a i 



V m1 U m1 r m V ml 



(60) 



a 



, ' m l 
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with the definition 

m— 1 

w ms = exp{i- (nks ~ 1)} 
1 k=i 

7T 



exp{i-(n fa - 1)} ( 61 ) 



7T — 

r m = exp{i- 22 ( n *T ~ X )} 

fc=m+l 



Under the transformation, the L-operator changes into |J: 

(jj)V m (62) 

where 

= Vm^Um} (g> V m \JJ m i 

U m ^ s = dia(w 
Substituting equation (|62|) into £#(//), we obtain 



where 



ex P{-^EK + ^-2)} = e-* 
4 i=i 



e 



- - e W {-*jj:(v]-T;)} = e-^ (65) 
4 i=i 

This is exactly equal to equation (46) with condition floTf ). This finished our proof. 

It is worthy to point out that Wadati et al also applied the Jordan- Wigner 
transformation to the L-operator and the R matrix. In the derivation of the hamiltonian 
of the Hubbard model, they impose the periodic boundary condition. But, they did 
not consider the relation between the boundaries. 

In the rest of this section, we give an another independent proof in terms of the 
hamiltonians. The periodic Hubbard model is 

L-l L 1 1 

H = - E (Ci/"'^ + a m,s a m+i,s) + (ji m] - -){n rni - -) 

m=l,s m=l,s (66) 

- ( a t.s a L,s + at,s a i,s) 



(63) 



e^TuOu) + e^T 22 (/i) + e ft T 33 (/i) + e^T u (ji) (64) 
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where, we have use the periodic condition 

a l+i,s — a t, s i a L+i,s = a>i, 8 - (67) 
Using the Jordan- Wigner transformation, we can obtain 

a m+l,T a «i,t + a m,T am +l,t — _ ( ff m+l ff m + Cr m ff m+l) 
a m+l,i a m,l + a m,l a m+l,l = ~( T m+l T m + T m T m+l) 



in L 1 

a i,] a L,] + att a i.T = ex P{y J2( a j ~ 2^ a L a t + a L°i ) 

in L 1 

a^a La + 4 a a u = exp{ — ^(r/ - -)}(r L r+ + r+rf) 

j'=i 



(68) 



The hamiltonian changes into 



L_1 ivr L 1 

rn=l Z j=l Z 

ivr L 1 

+ E (^ + i^ + W+) - exp{- Y>/ - -)}(r L - r+ + r+rf) (69) 

m=l Z j=l Z 



41 m=l 



2 r z 



Therefore, under the Jordan- Wigner transformation, the hamiltonian of the Hubbard 
model becomes one of the coupled XY model with twisted boundary condition 



ot» = exp{±^f:K-l)K d 
rt +1 = exp{±^ f>/ - l)}r= 

A 3=1 



j=1 (70) 



which is coincident with equations (57) and (58). 

7 Concluding remarks 

In this paper, we have found the eigenvalues of the Hubbard model and the coupled XY 
model with twisted boundary condition by using the analytic Bethe Ansatz method. 
We have shown how they are equal on the levels hamiltonian and the transfer matrix by 
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using the Jordan- Wigner transformation. The power expansion of log(t 9 (/i)) in terms 
of \x will give explicitly the infinite number of conserved quantities. In this level we 
claim the coupled XY model with twisted boundary condition is integrable. 

It is worthy to point out that we consider here only a set of special boundary 
conditions. It is not difficult to generalize to other kinds of twisted boundary conditions. 
For the open boundary, one should consider the solution of the reflection equations. 
This will be related to the surface critical behaviour of the system. 

Note: After finishing this paper, we was told by Prof. Wadati that Ramos and 



Martin |L4[ found the eigenvalue of the Hubbard model by using the algebraic Bethe 



Ansatz method, which recovers partly our results from the different approach. In 1 14 



they also notice the effect of boundary condition, but they did not discuss it in detail. 
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